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ABSTRACT 

RELATIVE PERFORMANCE OF DIFFERENT LMFBR LATTICE 
CONFIGURATIONS UNDER NATURAL CIRCULATION 

by 

John William Newton 



Submitted to the Department of Nuclear Engineering on 
August 24, 1970 , in partial fulfillment of the 
requirements for the degree of Master of Science 



Analytic comparisons of the performance of LMFBR open- 
hexagonal, triangular, and square lattice configurations 
under natural circulation flow conditions were made. 

Lattices having the same volumetric compositions were 
compared vjlth respect to reactivity change upon sodium 
voiding, pressure drop, and linear power output. 

A brief study of the relative advantage of the lattices 
from the standpoint of reactivity change upon sodium voiding 
of a coolant channel showed that the open-hexagonal lattice 
was preferable because of its larger coolant -channel radius. 

As regards hydraulic and thermal aspects, it was found 
that core pressure drops and linear power outputs were 
equivalent for the three lattices under botn natural 
circulation and forced convection, for both laminar and 
turbulent flow regimes. 

Consideration was given to developing more accurate 
relations for the hydraulic characteristics of the open- 
hexagonal lattice by modifying the traditional hydraulic 
diameter concept. Finally, some limitations of the open- 
hexagonal lattice as compared to the triangular lattice were 
discussed and necessary areas for future work were Identified. 
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Introduction 

1.1 Background 

In liquid-metal cooled fast-breeder reactors an area 
of particular concern is that of loss of heat removal 
capability and subsequent coolant boiling. The designs 
for most large LiDFBR's of the future include the capability 
for the removal of generated heat by natural circulation in 
the event of loss of pumping power. Most of these designs 
\iere based on the use of triangular pitch fuel pin lattices. 

1.2 Objectives 

The objective of this thc-sis is to determine the 
relative performance of different LMFBR lattice configura- 
tions under natural circulation with respect to the pov;er 
output which could be attained during normal operation and 
v;ith respect to safety considerations in the event of loss 
of primary pumping pov;er in a pumped circulation system. 
Figures 1-1 and 1-2 describe the lattice configurations 
under consideration, 

1.3 Discussion 

Advantages of natural, circulation operation Include: 

1. At full poxver operation primary pumps would 
be eliminated, a factor v;hich v;ould reduce cost 
of hardv/are and eliminate the concern of a loss 



of pumping power. 





Pig. 1.1 



Lattice configurations 



t11 




spacer 



a. Open-hexagonal cell 
unit cell 
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2. At full pov.'er operation a natural circulation 
system would, to some extent, be self -regulating; 
for a given channel, as power Increases/decreases, 
the core inlet and outlet temperature difference 
and the coolant mass flow rate would increase/ 
decrease. The parallel channel problem requires 
further investigation, however. 

3. In a pumped system natural circulation v;ould 
serve as a safety measure in that the loss of 
primary pumping pov/er v;ould not result in coolant 
boiling and burnout since a sufficiently large 
natural flow could carry away the heat quickly 
enough. 

Disadvantages of natural circulation operation Include; 

1. The plant will have to operate at a lovjer 
power output than would a plant with primary 
coolant pumps if design variables cannot 
adequately be optimized within various constraints 

2. The design of the system v/ould necessitate 
placing the intermediate heat exchanger far above 
the core in order to achieve the necessary height 

between the core and the heat exchanger for a 
desired thermal driving potential. In this 
position the intermediate heat exchanger could 
complicate access to the core from above and v/ould 
necessitate a larger (taller) outer containment 



structure 
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1.4 Choice of Design Constraints and Variables 

The loop-type system was chosen as the primary coolant 
system. Since the reactor system needs a high thermal 
driving head in order to attain the desired coolant flow 
rate, the intermediate heat exchanger needs to be above the 
core at a height which cannot practically be designed into 
a pot-type system. However, the pot-type system, through 
the thermal inertia of its large pool of sodium, provides 
ample margin for core cooling follovjing loss of pumping 
power . 

Sodium is by far the coolant of choice for recent LMPBR 
designs. Thus, in this study liquid sodium was chosen as 
the coolant. 

The type of spacer system chosen for the open-hexagonal 
lattice was the grid system. V/ire-wrap spacers would 
probably be more practical economically^ ^ ^ , but v/ith the 
center of the open-hexagonal lattice open, each rod would 
have support only from three wires in one complete spiral 
instead of support from six wires as in a triangular pitch. 
The grid spacer v;as chosen in order that the open-hexagonal . 
system would have the most and sufficient support against 
vibration. 

In the comparison of the three lattices several 
variables will be examined. Under natural circulation 
the linear povjer output depends on a combination of core 
temperature rise, fuel pin diameter and pitch, fuel pin 
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length, and the effective heat exchanger height to pz’lmary 
coolant system length ratio. Combinations of the variables 
v;ill be discussed with the use of typical values of current 
LMPBR fluid and mechanical parameters. 

1.5 Synopsis 

The comparison of the LMPBR lattices will first be 
examined from the standpoint of the additional leakage 
reactivity resulting from the heterogeneity of the core 
using the different lattices. An examination of the 
relative mass flow rates, core pressure drops, and linear 
power outputs vjill then be made follov.'ed by an example 
calculation of these quantities. In addition a more appli- 
cable expression for hydraulic diameter zvlll be sought for 
an open-hexagonal lattice. Finally, conclusions and recom- 
mendations v-;ill be proposed. 
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O’hapter 2 

Reactivity Effects of Voided Coolant Channels 

2.1 Reactivity Effect Components in a Sodium Voided Coolant 
Channel 

When a coolant channel is voided of sodium, a change in 
reactivity occurs. This reactivity change consists of three 
first order components: a change in the number of neutrons 

captured, a shift in the neutron spectrum, and leakage of 
neutrons. Detailed discussion of the total sodium-void 

( 2 ) 

effect may be found in publications by Hummel and Okrent' ’ 

f 3) 

and Monroe' 

Most calculations of the leakage component of the change 
in reactivity are based on the assumption that the reactor is 
homogeneous. In this chapter a rough calculation of the 
additional leakage reactivity effect due to the heterogeneity 
of the core v/ill be made. For a more detailed discussion of 
the leakage component of sodium-void reactivity in hetero- 
geneous cores, the reader is refered to a study made by 

( h.) ( 5 ) 

Lacapelle' ' v;hose method was based on the work of Benoist ' . 

The leakage component of the overall reactivity change will 
be compared in the open-hexagonal and triangular lattices 
showing that the larger hole of the open-hexagonal cell 
results in a greater negative reactivity than does the hole 
in the triangular cell. The calculations used in the com- 
parison are based on formulations derived in Behren's 
original paper "The Effect of Holes in a Reacting Material 
on the Passage of Neutrons" 



2,2 Leakage Reactivity Calculations for Open-hexagonal 
and Triangular Lattices 

Comparing Behren's Eq . (8) and Eq, (lO) for closely 
spaced, interacting holes in a medium having a long mean 
free path we have 



= JL 






'0 



A 



( 2 - 1 ) 



Inhere 



^ is the volume ratio of the holes to the material 
of the reactor, 

Q is a function of the shape of the hole, 

A is the mean free path in the solid, 

^hy ^ nuclear hydraulic radius of the cell, 

y Is the volume of the cell of material surrounding 

each hole, 



S is the surface area of each hole, 

2 

L^ is the mean square direct distance travelled by a 
neutron in the course of N collisions in a homo- 
geneous medium, 

2 

L is the mean square direct distance travelled by a 
neutron in the course of N collisions in a medium 
with holes. 

The change in reactivity is given by 



k 



= -(5^){k„ - l)/k„ . 

^Hom 



( 2 - 2 ) 



The definition of diffusion length gives 
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since 



= 



^a^fuel 



a^fuel 



,0 



V. 



tot 



(V 



fuel 



+ ''hole> 5^) 



(2-4) 



Equating Eqs. (2-1) and (2-3) 



D o 4 > 

(1 + f )— = (1 + 

Dq X 

or 

Qr, (b 

D - D„{1 + 4 ’) + . (2-5) 

^ 3(1 + <?) 

The change in the diffusion coefficient is 

6d = D - D„ = ?— -y --— - (2-6) 

Horn 2 

where 

^Hom homogeneous value of the diffusion 

coefficient 

Vhoie “ volume fraction of holes 

For a first approximation we XA)ill assume that the cells 
in our theoretical core are cylindrical. The volume of each 
cylindrical hole will be set equal to the actual volume of 
coolant channel associated v;ith a cell, and the volume of 
fuel material surrounding each hole will be set equal to the 
actual volume of fuel material associated with a cell. Refer 
to Fig. 1.2 for the unit cell configurations chosen. Note 
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that the cells are based on the coolant channel and not on 
a fuel rod as Is the more common convention. VJhen this 



is the same as the actual radius of the equivalent cylin- 
drical coolant cell. 

Refer to Fig. 1.1. a for a definition of parameters 
needed to find the hydraulic radius of the open-hexagonal 
cell; 



Refer to Fig. 1.1. b for the hydraulic radius of the 
triangular lattice: 




^hole ^tot ■^fuel * 




6J3 o 





(2-7) 



A 



hole 



A 



tot 



A 



fuel * 




( 2 - 8 ) 





(2-9) 



Inserting Eq . (2-6) In'co Eq. (2.2): 
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6k 

k 



= -( 



hole- 



3D. 



-)(■ 



nom 



■aO 




( 2 - 10 ) 



Using the following typical values for an open-hexagonal 
and a triangular lattice : 

Q = 4/3 (as tabulated in ref. (6)), 

'^hole = ^ “ '^fuel ^‘528 , 

p^ = 0.4l6 in., 

Ph = 0.340 in., 

d = 0.300 in. , 

^Hom ~ ~ 1*57 in. for a voided core, 

k ^ 1.5^ 

the ratio of reactivity changes is 



(6k/k)^ 

^ = 2 . 

(5k/k)^ 

Solving for the reactivity changes: 

(6k/k)^ = -0.0056 = -$1.70 

(5k/k)j^ = -0.0112 = -$3.40 

using p = 0.0033 2 $1.00. 

For a square lattice v;ith 

''hole = °-528 
d = 0.300 in. 

p = 0.388 in. 
s 




n 
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and 






1.42 



(6k/k)g^ = -$2.40 . 

From the above calculations we can see that the open- 
hexagonal lattice results in a higher negative reactivity 
from leakage in a voided coolant channel than do triangular 
or square lattices of the same volume fraction of fuel. 
Because of its higher negative reactivity the open-hexagonal 
lattice is preferable to both the triangular and the square 
lattices from the standpoint of keeping the effective 
increase in multiplication factor as lov; as possible upon 
sodium voiding of coolant channels. 

In the next chapter the relative mass flow rates, core 
pressure drops, and linear power outputs of the three lat- 
tices will be compared. These comparisons will be used in 
determining v;hlch lattice is preferable from hydraulic and 
pow'er removal aspects. 
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Chapter 3 

Expressions for and Comparison of Mass Flow Rates / Pressure 
Drops and Linear Power Output for Different LMFBR Lattice 

Configurations 

3.1 Discussion 

As was seen in Chapter 2 , the hexagonal lattice proved 
to be more advantageous than a triangular or square lattice 
from a leakage reactivity standpoint. In this chapter the 
three lattices will be compared with respect to some hydrau- 
lic and power aspects to see if the larger hole of the open- 
hexagonal lattice renders it "preferable to both the square 
and triangular lattices. In the present chapter the three 
lattices will be compared on an equivalent basis, namely, 
equal rod size and fuel volume fraction, with respect to 
pressure drops and linear power output. It will be assumed 
that all fluid properties are equal for the three lattices. 

As will be seen in the numerical examples of Chapter 4, 
the core pressure drops of the three lattice arrangements are 
essentially the same for equal fuel volume fraction and equal 
coolant velocities. The only differences in pressure drop 
are due to the differences in spacer arrangement of the 
different lattices, and these differences are small. There- 
fore, these small pressure drop differences can be neglected 
in a first order approximation of the relative linear power 
output for the three lattice arrangements. However, in the 



I 
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calculation of the absolute linear power output for any given 
lattice, the spacer pressure drop contribution to the total 
core pressure drop cannot be neglected since it constitutes 
as much as 45% of the total core pressure drop. 

3.2 Procedure for Calculation of Core and Spacer Pressure 

Drop 

The total pressure loss, AP^, through the reactor core, 
hence, through a unit cell flow path in the core, consists 
of three factors: pressure loss at bundle entrance and 

exit, APj,; frictional pressure loss without spacers, AP^^; 
pressure loss due to spacers, APg. It will be assumed that 
pressure loss at the bundle entrance and exit will be equal 
for all three lattices. This assumption is based on the fact 
that the entrance and exit openings are of equal total cross 
sectional area for the cores containing the different lattices 
and that the coolant velocities in the three cores are equal. 
By using equations similar to Eqs. (3-1) and (3-2), one can 
find that the entrance and exit losses for the different cores 
are the same. 

Since there is no experimental data on pressure drops 
through open-hexagonal rod bundles with and without spacer 
devices knovm to the author, the approach for calculating 

the semi-empirical pressure drop through an open-hexagonal 

(7) . . 

lattice will be as follows. Tong's method for finding 

the theoretical pressure drop through a channel with spacers 
is used to find theoretical values of APg for the triangular 



I 
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and open-hexagonal latti.ces. “Then a AP_ is found from 

O 

experimental data gathered and correlated by de Stordeur^^^ 
for triangular lattices. A ratio of experimental to theo- 
retical spacer pressure drops is then obtained for trian- 
gular lattices. The theoretical pressure drop for the open- 
hexagonal lattice calculated from Tong's formulae are then 
multiplied by the ratio found above to obtain an appropriate 
semi-empirical spacer pressure loss for the open-hexagonal 
lattice. This approach thus uses the existing experimental 
data of de Stordeur to normalize the calculated results- from 
Tong's method v^hich is applicable , unlike de Stordeur 's data 
to both triangular and open-hexagonal lattices. 



3.3 Theoretical Core and Spacer Pressure Drop 

According to Tong "the pressure losses in a rod bundle 
due to the rod spacers are the form drag type pressure loss. 
This type of pressure drop can be calculated by using the 
pressure loss coefficients at sudden contraction, and 

at sudden expansion, K^, as given by Kays and London 
Values of and are reproduced in Fig. 3.1. The pres- 
sure losses at sudden contractions and expansions are given 
as follows ; 

For a sudden contraction 



APc = 



2g^ 



(l-a ) + K 






per spacer. 



(3-1) 



1.0 



K 



0.8 



0.6 



0.4 



0.2 - 



0.0 



- 0.2 




Turbulent 



10,000 

5,000 

3,000 



Turbulent 

N = 3,000- 
5 , 000 ' 
10, 000 




-0.4 - 



0.0 0.2 0.4 0.6 

cr = Ag/A, 



0.8 



1.0 



Fig. 3.1. Pressure loss coefficients at sudden 

contraction , K , and at sudden expansion, K , as given 

c / q\ e 

by Kays and London''^-'. Taken from Ref. (7), Fig. 7. 
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For a sudden expansion 

2 2 

pV^ 2 

AP = - (1-a ) - K per spacer, (3-2) 

2g ®2g 

where 

a = A 2 /A^ is the ratio of the flow in the gridded region 
to the flow area away from the gridded region, 

V 2 is the velocity in the restricted cross-section 
region of the channel. 

Note that positive AP's are pressure losses. Refer to 
Fig. 3.2 for definition of flow areas. 

Therefore , 



AP^ - APg = APj^ — + (APg) , (3-3) 

cir 



where 



is the ratio of the friction factors in non- 



^cir circular to circular channels; see Fig. 3.3^^^^ 
Substituting an expression for the friction pressure drop 
in a duct into Eq. (3-3) gives 



4plv^ 



AP - AP = f ^ ( ) + AP + AP . (3-4) 

^ ^ 2g f . ® 

e ^c cir 



Substituting expressions for AP^ and AP^ from Eqs. (3-1) and 
(3-2) into Eq. (3-4) 



4piV^ f 






AP^ - a?e = fp 



D 2g f . 
e ^c cir 



( ) + (K^ + K^) 



n 



2g, 



(3-5) 




Region blocked 
by spacers 



Pig. 3o2o Flov; areas in a coolant cell, 
is the area avjay from the spacer region. A 
the free flovj area in the spacer region. 



cir 
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Fig. 3.3. Fully developed turbulent flow parallel to 
a bank of circular tubes of rods; points O based on friction 
coefficients calculated by Deissler and Taylor^ at Re = 10^; 
Reynolds number influence is small, and Nusselt number 
behavior is virtually the same as friction behavior. Taken 
from Ref. (lO), Fig. ( 9 -IO), which also shows that the ratio 
of Nusselt numbers is similar to the ratio of friction factors 
for non-metalic fluids. 
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where 



rig is the number of spacer grid sets in the rod bundle. 



Using 

fp = 0.046 Re for turbulent flow, 

Eq. (3-5) becomes 

0.092y°’^lp°* f 

AP^ - APp = D^*^g ^f . ^ 

e cir 



p (A3^/A2)^^1 

2g^ 



(K + K )n , 
e e g 



(3-6) 



where is found from V = co/pA^ (ft/hr) , and w is found 
from 

g3p^Z D^*^S^*^q 0.357 

0) = ( 2 — i) (Ib/hr) . (A- 8) 

0.092y'^*^Cp 

The quantity w is the average mass flow rate for natural 
circulation derived in Appendix A, n^ is the number of effec- 
tive rods per cell and q = n^q'l. 



3.4 Empirical Spacer Pressure Drop 

An empirical AP for a triangular lattice will now be 

( 8 ) 

found using de Stordeur's^ ' method. According to de 
Stordeur 



AP 



S 



2gAj_ 



(3-7) 



where 

C„ is the spacer drag coefficient 

O 

V is the velocity in the spacer region (ft/sec) 
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S is the spacer projected frontal area (ff^) 

2 

g is the gravitational acceleration (ft/sec ) 
is the flow area in the cell away from the 



gridded region (ft^) 



. 2 

p IS the density of the coolant (Ib/ft ) 

. 2 

AP„ IS the pressure loss due to one spacer (Ib/ft ) 



Using appropriate values of Vg, S, g, p, and a Cg 

taken from Fig. 3.4, an empirical spacer pressure drop, 
can be calculated from Eq. (3-7) for a triangular lattice with 
a honeycomb pin support grid. 

We can now obtain the ratio of the empirical spacer 
pressure drop to the theoretical spacer pressure drop for 
a triangular lattice with a honeycomb grid and set it equal 
to the ratio of the empirical to theoretical spacer pressure 
drop for an open-hexagonal lattice: 



An appropriate found from Eq. (3-8) for an 

open-hexagonal lattice after a open-hexagonal 

lattice has been calculated using Tong's method. Equation 
(3-6) then becomes 




(3-8) 



(AP^ - AP^)j^ 




) + 




(3-9) 



Drag Coefficient 
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Fig. 3.4, Drag coefficients for honeycomb grid spacer, 

(Dashed line is an extrapolated portion of the curve.) The 

quantity d is the spacer width. ' ' 
s 



- 31 - 



3.5 Mass Flow Rate and I, inear Power in Steady-State Thermal 
Circulation 

In steady-state thermal circulation the thermal buoy- 
ancy driving force, equals the friction pressure drop, 

APj^, for a pipe loop including the rod bundle section. See 
Appendix A for a derivation of mass flow rate and linear 
power relationships. For a rod bundle section the pressure 
losses through the spacers adds another component to the 
total pressure loss in the loop. The steady state pressure 
balance will then become 

f 

APg = APj^ + APg = fgAPj^ , (3-10) 

^cir 

where f^ is a multiplying ratio which accounts for spacer 
pressure drop. It is introduced in order to make the cal- 
culation of mass flow rates, pressure drops, and linear 

power output easier. The factor f is found to be approxi- 

s 

mately constant for the range of velocities investigated in 
any given lattice configuration. 

Substituting the relationships for AP^ from Eq. (A-3) 
and AP„ from Eq. (A-5) into Eq. (3-10) and solving for the 

X\ 

mass flow rate gives 

g$p^Z D^*^S^‘^n q'l 0.357 

w=( 7T9 ^ ' (3-11) 

0.092y^* Nc f 
«. p s 

where to is the mass flow rate for a unit cell. 

As shown in Eq. (A-3) the thermal buoyancy pressure gain 
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depends on the linear power and the mass flow rate, and the 
friction pressure drop depends on the mass flow rate. 

These dependencies can be seen in Fig. 3.5. In Fig. 3.5 
Eq. (3-11) is the locus of points along the AP_f line, and 
for any given the correspondting w can be taken from the 
graph. The dashed lines represent constant t^ - t^ lines. 
Since 



c At 




for any given core lattice geometry, any one of the three 
variables q' , w. At can be found, given the other two. 
Solutions for a natural circulation operating condition will 
be the point of intersection of the three lines representing 
AP for a given q', AP f , and At. Substituting Eq. (3-11) 
into Eq. (3-12) 

g3p^Z D^’^S^’^n 0.555 At 1.555 c 1.0 

g' = ( '■ ^ (— ) (3-13) 

0.092p-^' Nf n_ 1 

ft. s r 

Equation (3-13) is the relationship for the intersection of 

the AP_f line in Fig. 3.6 with a given At. The quantities 
1\ s 

(0 and AP„ are taken into account in Eq. (3-13). 

3.6 Comparison of Linear Power Output in Open-hexagonal, 
Triangular, and Square Lattices 

As shown in Fig. 1.2, the number of effective rods per 



cell for the three lattices is 



-33 




Fig. 3.5c Illustrative graph of solutions for natural 
circulation operating conditions. The intersection of the 
three lines defines an operating, condition. The AP^^f^ line 
is defined by Eq . (A-6). lines are defined .by Eq. (A-3) 

The At lines are governed by Eq . (A-11). 
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open-hexagonal - n^ = 2 



- = 1/2 



- n = 1 
r 



triangular 
square 

For this comparison it will be assumed for a first order 
approximation that the fluid properties of the coolant, 

At, Z , f , and N are equal for the three lattices. The 
three lattices will be compared by taking a ratio of the 
linear power outputs of the three lattices based on equal 
fuel volume fraction and equal pin diameters. The variables 
will be the equivalent diameter, D^, of a flow cell, the 
cross-sectional flow area, S, of a cell, and the number of 
effective fuel rods per cell, n^. As shown in Appendix B, 
for equal the hydraulic diameters of the three lattices 
are equal. By definition 



D = 



4 X cross sectional area 
wetted perimeter 



4S 



w 



or 



S = 



P D 
w e 



Therefore , 



S, = 






udD 
€ 

2 

TfdD 
€ 

8 

udD 



sq 



4 
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Using Eq. (3-13) for q ' , 



q' 



g3p^Z 0.555 At 1.555 c 1.0 

^ — "0 2 ' (3-13) 

0.092^‘^Nf n 1 



we obtain the following ratios of linear power output: 





^n ),‘^*^^^(l/n ), 
r h ' r h 

8 ,0.555/,, , 

n ) (1/n ), 

r ' ' r h 



1.555 

1.555 



( (TTdD^/2)^*® (2) )°*^^^ (1/2)^*^^^ 

( (TTdD^/8)^*® (1/2) )°*^^^ (1/0.5)^*^^^ 

( (TTdD^/2) (2) )°*^^^ (1/2)^*^^^ 

q' ((TTdD /4)^*®(1))^-^^^(1/1)^-^^^ 

o QJ w 




It is seen that the linear power output of the three lat- 
tices under natural circulation is the same for equal fuel 
volume fraction and equal pin diameters. This is a direct 
consequence of the fact that the hydraulic diameters are the 
same. It is also easily shown that power removal capabilities 
under forced convection laminar or turbulent flow conditions 
are also the same for the three lattices (see Appendix D) . 

In Chapter 4 some typical values of fluid and lattice 
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parameters will be used to calculate the pressure drops and 
mass flow rates for a range of linear power output. Using 
typical core and primary system parameter values, a q' for 
natural circulation will be calculated to see how this value 
compares to current linear power ratings for forced convection 
cooling. 
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Chapter 4 

Example Calculations of Mass Flow Rate, Pressure 
Drop and Linear Pov;e r Output' ' 

4.1 Discussion 

In Chapter 3 the relationships for the mass flow rate, 
the pressure drop, and the linear power output per cell 
were derived. In this chapter typical values for lattice 
parameters will be used to solve for the mass flov; rate, 
the pressure drop, and linear power output for an open- 
hexagonal lattice and for a triangular lattice of equal fuel 
volume fraction, V^, and equal pin diameter. As was seen 
in Chapter 3^ the Ij.near pov/er output for triangular, open- 
hexagonal and square lattice arrangements are equal for 
equal fuel volume fraction, for equal core temperature drop, 
and for equal fluid properties in the three lattices. This 
equality is based on the condition that the pressure drop 
through the flo;i/ cells of the three lattices is equal since 
coolant velocities are equal in the three lattices. This 
condition is not strictly accurate because of the differ- 
ences in the spacer arrangements. These differences intro- 
duce only small second order errors and can be neglected 
in a first order calculation of the relative linear pov/er 
output. However, since they cannot be neglected in more 
precise calculations of the mass flow rate and linear power 
output of any given lattice arrangement, we will include 
calculations of these effects and tables of values of mass 



I 
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flow rates, velocities and core pressure drop for a range 
of linear pov;er output from one to tv;enty kilowatts per 
foot for the different lattices. 

4.2 Calculation of Flov/ Parameters and Theoretical Spacer 
Pressure Drop in the Open-hexagonal Plow Cell 

Refer to Figs. 1.1. a and 1.2. a for definition of the 
parameters required to calculate the expression for the 
of an open-hexagonal lattice. The values of the para- 
meters used to calculate (Vl), are d = 0.300 in., 6. = 0.040 

' I ' h p h 

in., and p^ = 0.34o in., from which 

■Jry[3 d p 

(vph = — 

9 

The value for d is about as large a value as is being 

P 

considered for present LMFBR's fueled with uranium dioxide. 
The value for 6^^ is about as small a pin to pin spacing as 
is being considered. 

The following values will be used in calculating mass 

flow rates, pressure drop, and linear power output for an 

open-hexagonal lattice flow cell; 

n^ = number of pins per cell = 2, 

A t= cross sectional area of a fuel Din = 7rd^/4 
P , P * P 

== 0.00049 ft , 

d = assumed spacer web thickness = 0.00125 ft, 
s 

n = number of spacer segments per cell = 9 
s 

Ig = length of one spacer segment = p^/3/3 = 0.0164 ft. 
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Ag = cross sectional 
= 0.000184 ft^. 



area of spacers per cell = n 1 d 

s s s 




A 



1 



flow area away from the gridded region = A^ 

0.0011 ft^. 



n A 
r p 



Ao = flow area in the gridded region = A, - A 

^ r> 2 Is 

= 0.000918 ft (see Fig. 3.2 for definition of flow 
regions) , 

^ = Ag/A^ = 0.834. 

Equation (3-11) can now be used to calculate the variation 
of flow rate with power output. 



cv= (. 



g3p 



o'l 



e e 



0.092m-^*^N c f 
e p s 



r- :^0.357 



(3-11) 



with the following additional constants: 
g = 4.16 X 10^ ft/hr^, 

P = 1.6 X 10"^ 

P =50 Ib/ft^, 

1^ =0.5 Ib/ft-hr, 

c =0.3 BTU/lb-°F, 

(Ze/^e) = 0.5> 

= ( 3 >Gp^ - Trd^)/TTdp'= 0.0281 ft, 

S = = 0.0011 ft-, 

=2 rods, 

1 = 4 ft, 

=1 (see Eq. (3-10) for a definition of f ) . 

s 





I 




-4o- 



Tlius, for example, for q' - 10 kv//ft = 34130 BTU/trr-ft, the 
average mass flow rate for an open-hexagonal lattice flow 
cell Is = 602 Ib/hr. Table IV. 1 lists values of mass flow 
rate, coolant velocity in the ungridded region of a flow cell 
and corresponding Reynolds number for the open-hexagonal 
lattice. The coolant flovj is seen to be turbulent for the 
range of linear power output considered. 

The following computations are sample calculations of 
the average coolant velocities in and away from the gridded 
region of an open-hexagonal flow cell for the q* = 10 kw/ft 
case. These velocities will be used in finding the pressure 
drop due to the fuel pins and the spacers in the open-hex- 
agonal flow cell. The quantity is the average coolant 
velocity away from the gridded region: 

(Vi)h = ^/p(Ai)h = 10,966 ft/hr . 

The quantity Vg is the average coolant velocity in the 
4 gridded region: 

(^2^h = = 13,164 ft/hr . 

The quantity Re^ is the Reynolds number for the coolant 
away from the gridded region: 

(Rei)h = "" 30,815 . 

The quantity is the Reynolds number for the coolant 
in the gridded region: 



I 
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Table IV . 1 

Flow Parameters for the Open-hexagonal Lattice 



q' (kw/ft) 


(ib/hr) 


(ft/hr) 


Re^ 


1 


265 


4,820 


13,544 


2 


340 


6,173 


17,347 


3 


392 


7,135 


20, 049 


4 


435 


7,906 


22,217 


5 


471 


8,562 


24,060 


6 


503 


9,138 


25,678 


7 


531 


9,655 


27,130 


8 


557 


10,126 


28,455 


9 


581 


10,561 


29,677 


10 


6o4 


10,966 


30,815 


11 


624 


11,346 


31,881 


12 


644 


11,704 


32,887 


13 


662 


12,043 


33,840 


14 


680 


12,366 


34,748 


15 


697 


12,674 


35,614 


16 


713 


12,969 


36,444 


17 


729 


13,253 


37,242 


18 


744 


' 13,526 


38,009 


19 


758 


13,790 


38,750 


20 


772 


1^,045 


39,466 
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- Dgp(V2)j^/u. - 36,990 

From Fig. 3.3 for (T = 0.834, K = 0.14 and K = -0.02. 

c e 

Using Eqs. (3-3) and (3-6), the value for the total flow 
cell pressure drop is calculated: 



(Ap^ - 



(3-3) 



cir 



0 . 2 , _ 0 . 8 ,, 1.8 



0.092n^*^lp"‘ X*" f, 

(APj - APj.) = 



So 



cir 



P(A^/A2)^V^ 



(K K )n (3-6) 
2g. e c g 



(AP^ - ^Pg^h " ^ 10"S]^*^(1.02) + 0.052 X 10“°vJ 



^-6„2 



(AP^ - APg)j^ - 30.7 Ib/ft^ . 



Table IV. 2 lists values of AP , (APp,^; )j^, ( (AP„)^j^) 

(Ap^ - APg)p^ for a range of linear power output. Since 
Reynolds numbers and velocities are equal for triangular and 
open-hexagonal lattices, Ap ' s v/ill be equal for the tv.'o 

n 

lattices. 



4.3 Calculations of Flov; Parameters and Theoretical Spacer 
Pressure Drop in the Triangular Flow Cell 

Refer to Figs. 1.1. b and 1.2.b. for definition of the 

parameters required to calculate the expression for the 





I 
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Table IV. 2 



Pressure Drop in the 



Open-hexagonal Lattice Flovj 



Cell 



[ 

1 (ft/hr) 


AP^ ' 


(Ih/ft'^) 


((^^sUh^h 

(Ih/ft^) 


(Ib/ft^) 1 


1 4,820 


5.45 


5.56 


1.21 


6.77 1 


1 6,173 


8.51 


8.68 


1.98 


10.67 ] 


7.135 


11.05 


11.27 


2.65 


13.90 1 


1 7>906 


13.29 


13.56 


3.25 


16.80 1 


1 3^ 


15.34 


15.65 


3.81 


19.46 1 


r 9^138 


17.25 1 


17.59 


4.34 


21.94 1 


1 9,655 


19.04 


19.43 


4.85 


24.27 J 


1 1 o . T p6 


20.75 


21.16 


5.33 


26.50 1 


o < 

• 'i 

1 

H ' 

— 


22.38 


22.83 


5.80 


28.63 1 


1 10,966 


23.95 


24.43 


6.25 


' 30.68 ] 


31,346 


25.^6 


25.97 


6.69 


32.67 1 


11,704 


26.93 


27.47 


7.12 


34.59 1 


1 12.043 


28.35 


28.92 


7.54 


36.46 


1 12,366 


29.73 


30.33 


7.95 


38.28 1 


12,674 


31.08 


31.70 


8.35 


40.05 1 


12,969 


32.39 


33.04 


8.75 


41.79 1 


1 13.253 


33.68 


34.35 


9.13 


43.49 


1 13,526 


34.94 


35.64 


9.51 


45.15 


1 13,790 


36.18 


36.90 


9.89 


46.79 


1 ]L 4 , 045 


37.39 


38.14 


10.26 


48.39 
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of a triangular lattice. Using the same fuel volume fraction 



and fuel pin diameter as v;ere used for the open-hexagonal 
lattice, 

(vpt = -~(^)^ = 0-^72 . (4-2) 

2^^3 p^ 

Solving for the pitch and the pin to pin spacing gives: 
p^ = 0,4l6 in. 

6^ = 0.116 in. 

The follov;ing values will be used in calculating mass flow 
rates, pressure drop, and linear pov/er output for a trian- 
gular lattice flow cell: 
n^ = 1/2, 

Ap = 0.00049 ft^, 
d = 0.00125 ft, 

S 

"s = 1-5' 

1 = 0.02 ft, 

s 

A^ = 0.0000375 ft^, 

Arj, = nT 3 p ^/4 = 0.00052 ft^, 

A^ = 0.000275 ft^, 

Ap = 0.000237 ft^ (see Fig. 3.2 for definition of 
flow regions), 

C" = Ap/A^ = 0.865. 

Equation (3-11) can again be used together v/lth the follov/ing 
additional constants: 

g = 4.16 X 10^ ft/hr^, 

p = 1.6 X 10“^^ °f"^. 
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p =50 Ib/fb^, 

li, =0,5 Ib/ft-hr, 

c =0.3 BTU/lb-°F, 

(Z/Ng) = 0.5, 

Dg , = (2>/3p^ - •n-dp)/7rdp = 0.0281 ft, 

S = = 0.000274 ft^, 

n = 1/2, 

1 = 4 ft. 



Thus, for example, for q’ =10 kw/ft = 34130 BTU/hr-ft, the 
average mass flovj rate for a triangular lattice flow cell 
is a; = 150.5 Ib/hr. 

Table IV. 3 lists values of mass flovj rate, coolant 
velocity away from the gridded region, and the corresponding 
Reynolds number for the triangular lattice. Note that the 
corresponding coolant velocities and Reynolds numbers are 
the same for the open-hexagonal lattice since the fuel 
volume fractions are equal. Note also that the mass flow 
rates for the triangular lattice flov/ cell are one fourth 
of the corresponding mass flow rates for an open-hexagonal 
lattice flow cell. The differences are due to the differ- 
ences in flov7 area, A^, per cell and in the number of 
effective fuel rods, n^, per cell. The mass flow rates 
per unit flow area, equal for the tv;o lattices. 

The following computations are sample calculations of 
the average coolant velocities in and away from the gridded 
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Table IV. 3 

Flo\-j Parameters for the Triangular Lattice 



q * (kw/ft) 


uo (ib/hr) 


(ft/hr) 


Re^ 


1 


66 


4,820 


13,5^^+ 


2 


85 


6,173 


17,347 


■ 3 


98 


7,135 


20,049 


4 


108 


7,906 


22,217 


5 


Il8 


8,562 


24,060 


6 


125 


9,138 


25,678 


7 


132 


9,655 


27,130 


8 


139 


10,126 


28,455 


9 


145 


10,561 


29,677 


10 


150 


10,966 


30,815 


11 


155 


11,346 


31,881 


12 


161 


11,704 


32,887 


13 


165 


12,043 


33,840 


14 


170 


12,366 


34,748 


15 


174 


12,674 


35,614 


16 


178 


12,969 


36,444 


17 


182 


13,253 


37,242 


18 


186 


' 13,526 


38,009 


19 


189 


13,790 


38,750 


20 


193 


14,045 


39,466 



-Ml 

region of a triangular flov/ cell for the q ' = iO ku'/ft case. 
These velocities will he used in finding the pressure drop 
due to the fuel pins and the spacers in the triangular 
lattice flow cell. The procedure is the same as that used 
in finding the pressure drop in the open-hexagonal lattice 
flow cell in Section 4 . 2 . 

' (V^)^ = ^/p(A^)^ = 10,966 ft/hr, 

(Vg)^ = 12,701 ft/hr. 

Note that since the fuel volume fractions are 

equal, hut is slightly less than since the grid 

configurations are different for the two lattices. 

(Re^^)^ = = 30,815 

(Re2)t = = 35,688 

From Fig. 3.3 for = O.865, K = 0.12 and K = - 0 . 03 . 

c o 

Again using Eqs. ( 3 - 3 ) and (3-6), the value for the total 
flow cell pressure drop is calculated: 

(AP^ - APg)^ = 1.28 X 10“^^ *^(1.095) -t 0.036 X 10“^V^ 

(APt - APg)^ = 30.54 Ih/ft^ . 

The lattice pressure drop (AP^ - APj-)^^ and (AP^ - AP^)^ 
are very close in value. This is due to the offsetting 
effects of the values of the factors and the values 
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of the spacer pressure drop on the values of - ^Pg) 

of the tv;o lattices. Table IV, 4 lists values of Ap , 

^ range of 

clr 

linear power outputs. 



4.4 Empirical Spacer Pressure Drop for the Triangular and 
Open-hexagonal Lattices 

The method of de Stordeur for finding the empirical 
spacer pressure drop In a triangular lattice was presented 
In Section 3.4, Inserting the follov.’lng values Into Eq, 
(3-7)j here repeated, and solving for ((AP^)^^) 4 ., 

S GTTl X# 



AP 



s 




(3-7) 



g = 4.16 X 10^ ft/hr^, 

= 0.000274 ft^, 

= 12,687 ft/hr, 

S = A = 0.0000375 ft^, 

S 

p =50 Ib/ft^, 

C =2.33 (see Fig. 3,4 vjhere Re = d V p/p, = 1,586), 
s s s s 

gives 

(^^^s^em^t = IL/ft^-grld .. 



Using the theoretical spacer pressure drop calculated for 
a triangular lattice In Section 4,3: 






th't 



). = 



4.34 Ib/ft^ 
5 grids 



= 0,868 Ib/ft^-grld 
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Table IV, 4 

Pressure Drop in the Triangular Lattice Flov; Cell 



^1 


[> 








(ft /hr) 


(Ib/ft^) 


(Ib/ft^) 


(Ib/ft^) 


(Ib/ft®) 


4,820 


5,45 


5.97 


0.84 


6,80 


6,173 


8.51 


9.32 


1.37 


10.69 


7,135 


11.05 


12.10 


1.84 


13.93 


7,906 


13.29 


14.55 


2.25 


16.80 


8,562 


15.34 


16,80 


2.64 


19.44 


9,138 


17.25 


18.89 


3.00 


21.90 


9,655 


19.04 


20.85 


3.36 


24.20 


10,126 


20.75 


22.72 


3.69 


26.41 


10, 561 


22.38 


24.51 


4.02 


28.52 


10,966 


23.95 


26.23 


4.34 


30.57 


11,346 


25.46 


27.88 


4.64 


^ 

32.51 


11,704 


26.93 


29.49 


4.93 


34.42 


12,043 


28.35 


31.04 


5.22 


36.26 


12,366 


29.73 


32.55 


5.50 


38,06 


12,674 


31.08 


34.03 


5.78 


39.81 


12,969 


32.39 


35.47 


6.06 


41.53 


13,253 


33.68 


36,88 


6.32 


43.20 


13,526 


34.94 


38.26 


6.58 


44.84 


13,790 


36.18 


39.62 


6.84 


46.46 


14,045 ' 


37.39 


40.94 


7.10 


48.04 



Therefore j 



(Ap ) 
s "t h 



The pressure drop for the open-hexagonal lattice is 



o.o92n°'hp°’®vh® f 

- '^PE^h r?— — ^(— ) 

So ^clr 



(K. -I- K„)n„ , 



+ ( 7 -% ■■). 



2g, 



e c' g 



(AP^ - APg)j^ = 24.4 + 3 . 59 ( 6 . 25 ) = 46.8 Ib/ft . 



Calculating the factor, f , gives (see Eq. (3-10)) 

s 



AP„f/f . + (Ap ) 

(f^)^ = — ^ ^ = 1.96 . 



AP. 



R 



Table iV.5 lists values of V , Re , C , ( (AP )^„)+., 

s s s s Gin u 

finding (f ), for a range of linear power output from one 

S Al 

to twenty kilowatts per foot. (Ap^)^ and (^^^^/^cir^h 
found in Table IV. 2. The factor (fg);^ is relatively con- 
stant, v;ith an average value of 1.95. The quantity 
(AP^ - ^^E^h "" ^^Vs^h tabulated in Table IV. 5. 

For more accurate pressure drop calculations, (Tg)j^ 
should be used in Eq . (A-8) fo.r finding the mass flov; rate 
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Table IV. 5 

Flow Parameters and Pressure Drop in the Spacer Region' 
for the Triangular Lattice Flov; Cell 



q' 


V 


Re 


c 


((^P ) ) 4 - 

' s ern't 




S 


s 


s 


(kv;/ft ) 


(ft/hr) 


— 


— 


Ib/ft^-grid 


1 


5,577 


697 


2.68 


0.68 


2 


7,142 


893 


2.56 


1.07 


3 


8,255 


1,032 


2.51 


1.41 


4 


9,147 


1,143 


2.47 


1.70 


5 


9,906 


1,238 


2.44 


1.97 


6 


10,573 


1,322 


2.41 


2.21 


7 


11,171 


1,396 


2.39 


2.45 


8 


11,716 


1,464 


2.38 


2.69 


9 


12,219 


1,527 


2.36 


2.90 


10 


12,687 


1,586 


2.35 


• 3.11 


11 


13,127 


1,641 


2.34 


3.31 


12 


13,542 


1,693 


2.33 


3.51 


13 


13,934 


1,742 


2.32 


3.70 


14 


14,307 


1,788 


2.31 


3.89 


15 


14,664 


1,833 


2.30 


4.07 


16 


15,005 


1,876 


2.29 


4.24 


17 


15,334 


1,917 


2.28 


1 — 1 

• 


18 


15,650 


1,956 


2.27 


4.57 


19 


15,955 


1,994 


2.26 


4.73 


20 


16,250 


2,031 


2.25 


4.91 
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Table IV . 5 continued 



(Ib/ft®) 




Ib/ft^-grid 




(Ib/ft®) 


0.84 


4.10 


4.96 


1.93 


10.63 


1.37 


3.92 


7.76 


1.93 


16.59 


1.84 


3.83 


10.15 


1.94 


21.55 ' 


2.25 


3.78 


12.27 


1.94 


25.92 


2.64 


3.73 


14.20 


1.95 


29.90 


3.00 


3.68 


15.99 


1.95 


33.64 


3.36 


3.65 


17.72 


1.95 


37.13 


3.69 


3.64 


19.39 


1.95 


40.46 


4.02 


3.61 


20.92 


1.95 


43.64 


4.34 


3.59 


22.44 


1.96 


46.70 


4.64 


3.58 


23.92 


1.96 


49.65 


4.93 


3.56 


25.36 


1.96 


52.51 


5.22 


3.55 


26.74 


1.96 


55.28 


5.50 


3.53 


28.07 


1.96 


57.97 


5.78 


3.51 


29.34 


1.96 


60.60 


6.06 


3.59 


30.62 


1.97 


63.16 


6.32 


3.48 


31.80 


1.96 


65.68 


6.58 


3.47 


33.00 


1.96 


68.13 


6.84 


3.45 


34.16 


1.96 


70.55 


7.10 


3.45 


35.44 


1.97 


72.91 
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of the open-hexagonal lattice flov/ cell. Coolant velocities 

and cell and spacer pressure drop vjould be recomputed to' 

find a corrected This process v;ould be continued 

until succeeding vjere essentially constant. 

The factor, f , for the triangular lattice is calcu- 
s 

lated by the same procedure as was used to calculate 

The average value of (f^)^is l.'Jk, The factor, (f^)^, is 

not equal to (f ), since the spacer pressure drop in the 

S -H 

tv;o lattices are different, a condition due to the differ- 
ence in spacer configuration in the two lattices. 



4.5 Linear Povjer Output for the Open-hexagonal Lattice 
Using Eq. (3-13), 



(. S Pp""ePV y 'V)0.555(i^)1.555(!a)l-0 



0 . 092 (i‘^*^N f 
G S 



(3-13) 



n 



the average linear power for the open-hexagonal lattice 
is calculated using the following constants: 
g = 4.16 X 10® ft/hr^, 

■ p = 1.6 X 10“'^ °f"^, 

p =50 Ib/ft®, 

p, =0.5 Ib/ft-hr, 

Cp =0.3 BTU/lb/’?, 

1 kw = 3413 BTU/hr . 



In terms of several system variables, the linear pov/er 
output can be expressed as 
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22.4({^)0-555{^)0.555(,fi,1.55VSeS(J_)1.0) (^.3, 

N_ f_ 100 ® n 1 ' 

r 



where At is In ^F, is in inches, S is in square inches, 
and 1 is in feet. Inserting the following typical values 
into Eq. (4-3); 

Z /N =0,5 (i. e,, the ratio of the source/sink 

elevation difference and the effective primary 
circuit length). 

At = 400 ^P, 

2 

S = 0,1587 in, (the cross sectional flov/ area for 
the open-hexagonal flow cell in these calcu- 
lations), 

^ = 1 - 95 . 

Dg = 0.337 in. (the equivalent diameter for the 

open-hexagonal flow cell in these calculations), 

n^^ = 2 rods, 

1 = 4 ft, 

the average linear pov?er output for the open-hexagonal 
lattice (and for the triangular and square lattices) is 



= 0.873 kw/ft . 

The above value for q' is low compared to the required full 
power output of a practical reactor, vjhlch averages around 
10 kw/ft. On the other hand, 8 . 73 / of full power is a 
respectable rating in terms of decay heat removal. Sugges- 
tions for increasing the value of the linear power output 
achievable v;ith natural circulation flow are made in 
Chapter 6. 
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The mass flo;v rates, o>, coolant velocities, V, and 
friction pressure drop, calculated in the preceding 

example were based on an = 1, Table IV, 6 lists the mass 
flow rates per unit flov; area, total friction 

pressure drop, Ap^^f^, calculated using = 1.95 and 

(f ), = 1.7^ in Eq, (3-11), The tube bundle friction 
pressure drop, Ap^^, was calculated from 



APj^ = 1.28 X 10"^V^*® Ib/ft^ , 



where 



to 




Figure 4.1 shovjs the natural circulation operating 
conditions for the open-hexagonal and triangular lattices. 
Using a At = 400 °F, the linear pov;er output of the two lat- 
tices vary by the factor (l/f Since = 

4(A^eii)t ^h ^ horizontal axis of Fig. 4.1 is 

in units of ^/^qqii ease of comparison of the perfor- 

mance of the two lattices. Figure 4.1, v;hich is based on 
the values in Table IV. 6, shows that for equal total friction 
pressure drop and equal At, the triangular lattice has a 
slightly higher linear power output than does the open- 
hexagonal lattice. 
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Table IV. 6 

Mass Flow Rates per Unit Plow Area and Total Friction 
Pressure Drop as Functions of a Calculated f in the Open- 
hexagonal and Triangular Plov; Cells 



q' 

(kw/ft) 


^^^^^^cell^h 

(Ib/hr-ft^) 


(''Vs'h 

(Ib/ft®) 


^‘^'^^cell^t 

(Ib/hr-ft^) 


(Ib/ft^) 


1 


190, 000 


6.9 


198,000 


6.6 


2 


243,500 


10.8 


253,500 


10.4 


3 


281,000 


14.0 


292,500 


13.5 


4 


311,500 


16.8 


324,500 


16.2 


5 


337,500 


19.4 


351,500 


18.7 


6 


360,500 


21.8 


375 , 000 


21.1 


7 


380,500 


24.1 


396,000 


23.2 


8 


399,000 


26.3 


415,500 


25.3 


9 


416,000 


28.4 


433,500 


27.3 


10 • 


432,500 


30.4 


450,500 


29.2 



in (kw/ft) 
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Fig, 4,1, Natural circulation operating conditiona denoted ty O for the open- 
hexagonal and triangular lattices. 
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Chapter 5 
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Equivalent Diameter for the Open-hexagonal Lattice 

5.1 Definition of Equivalent Diameter 

Whenever flow in noncircular ducts is considered, it 
is necessary to employ in the calculation of the Reynolds 
number and friction factor some characteristic dimension. 

In circular ducts the characteristic dimension is the diameter, 
D. In noncircular ducts the concept of equivalent diameter, 

D^ , (also referred to as hydraulic diameter) applies when 
the duct does not appreciably vary from circular. The 
equivalent diameter is defined as 

D = 4 ^ cross sectional flow area _ 4A (5-1) 

^ wetted perimeter 

The equivalent diameter is the characteristic length which 
would give the same ratio of fluid passage volume to wetted 
area for circular and noncircular ducts. 

5.2 Conditions for a Laminar Flow Region in the Open- 
hexagonal Flow Cell 

In a noncircular duct the concept of equivalent dia- 
meter applies for turbulent or laminar flow. However, the 
values for the Reynolds number and friction factor are 
different for laminar and turbulent flow. In a unit open- 
hexagonal cell (refer to Fig. 1.2. a) the flow passage is 
approximately circular; therefore, the concept of equivalent 
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diameter should apply. In Chapter 4 the Reynolds number and 
friction factor for the open-hexagonal cell were calculated 
assuming fully turbulent flov/ in all parts of the cell. 

These numbers were computed using the average cell rate of 
flow which was turbulent. For turbulent flow 

fp = 0.046Ri®*^ . 

However, it is possible that the flow is laminar in some 
sections of the flow cell, namely, in the narrow flow area 
between two adjacent fuel pins. If this is the case, then 
the Reynolds number and friction factor are different in 
the laminar flow area from the Reynolds number and friction 
factor in the rest of the flow area. For laminar flow 

fp = 16Re^ . 

In order to determine if the flow in the region between 
adjacent fuel pins is fully laminar, a rough calculation of 
the laminar sub-layer thickness around the fuel pins will 
be made. The pin to pin spacing is 0.040 inches; therefore, 
if the laminar sub-layer around each pin is as much as 0.020 
inches, the. region may be considered to have laminar flow. 

If this be the case, the open-hexagonal flow cell can be 
divided into two channels, one having turbulent flow and 

one having laminar flow. This division of the flow cell 
assumes no inter-channel mixing and little or no heat trans- 
fer from the laminar region to the turbulent region. This 
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assumption is not good if there is a significant degree of 
inter-channel coolant mixing, which in fact probably occurs 
in the present case since the channels are not isolated. 

This assumption breaks down also if the heat transfer from 
channel to channel is good. 

If it is found that open-hexagonal cells can be divided 
into two channels, the coolant flow rates would have to be 
recomputed using two flow regions per cell in order that 
the power output could be recomputed. Since the pressure 
drop through the subassembly, will be the same for 

all flow regions, the coolant velocity in the laminar channel 
will be less than the coolant velocity in the turbulent 
channel in accordance with: 



(4f 



(AP 

LpV 

D 2 
e 



SA^ Lam 
2 

Lam 



(AP 

(4f 



SA^ Tur 

hpyl) 

D 2 ' Tur 
e 



r 



From this last relation it follows that 



V. 



Lam 



< V, 



Tur 



(5-2). 



since the friction factor for fully laminar flow is greater 
than the friction factor for fully tu.rbulent flow in a 
smooth pipe or duct, and the laminar region is approxi- 
mately less than or equal to the turbulent region D^. 
Calculations of the equivalent diameters for the two regions 

showed (D )., ^ (D in the .open-hexagonal flow cell, 

e Lam e Tur 



Tur 
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''hom ' ''lur' “Lam ^ <»Tur' “Lam 

flow rate in the laminar channel of the flow cell, and (0^ 

is the mass flow rate in the turbulent channel of the flow cell. 

If it is desired to maintain the core At at a certain amount, 

then (cj^. „ + _) < w where io is the flow rate based on 

Lam Tur sr sr 

a single region of fully turbulent flow in the entire flow 

cell. The sum of w, and is less than w because 

Lam Tur sr 

10 ,. „ in the laminar channel is less than . Since c and 

Lam Tur p 

At are constant and ) < w , then . will have 

Lam Tur sr' ^Tot 

to be decreased from the value calculated in a single region 
(turbulent) flow correlation (^ipQ^^sr ~ degrad- 

ing the power output of the entire core. 



5.3 Calculation of Laminar Sub-layer Thickness in Turbulent 
Flow 

The thickness of the laminar sub- layer around the fuel 

pins will be calculated to a first order approximation using 

( 12 ) 

equations given by Schlichting . The laminar momentum 
sub-layer thickness, 6^, for turbulent flow in pipes is 
given by Schlichting ' s Eq. (20.15a) as 




(5-3) 



The friction velocity, v^, is given by his Eq. (20.8) as 

V. = (T ^ ,0.03955iT-V^VVl/4,l/2 _ , 5 .,) 



T is the shear stress at the wall, - 
o 



where 
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u is the average velocity in the pipe, 

V is the kinematic fiscosity of the fluid, 
d is the pipe diameter. 

Using the following values from the example calculations 

in Chapter 4: 

d = (D ), = 0.0281 ft 

e h 

u = = 10,966 ft/hr 

V = y/p = 0.01 ft^/hr 

the friction velocity is v* = 598 ft/hr, and the laminar 
momentum sub-layer thickness ~ 0.001 in. The thermal 
boundary layer thickness on the other hand is ~ k/h. 

Using the following approximate values for liquid sodium: 
k = 40 BTU/hr-ft-^F, 
h = 20,000 BTU/hr-ft^-®F, 
the thermal boundary layer is ~ 0.002 in. 

These values of 6^ and 6^ are considerably less than 
0.020 inches, the value necessary to establish a laminar 
flow region between adjacent fuel pins. Therefore, the flow 
in all regions of the open-hexagonal flow cell can be 
considered to have turbulent flow. Consequently, the open- 
hexagonal flow cell can be considered to have only one flov/ 
region, and an equivalent diameter, defined by Eq. (5-1), 
can be calculated for the entire cell. 

The above calculations of 6^ and were only approxi- 

mate. A more accurate computer calculation using a code 



can be made 



(13) 

such as General Electric's Velvet-II Code 
of the velocity profiles around the fuel pins. However , 
there is no reason to expect that more sophisticated cal- 
culations would alter the above conclusions. 
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C'hapter 6 

Summary and Conclusions 

6 . 1 Summa ry 

The work reported in the preceding chapters has compared 
the open-hexagonal, triangular, and square lattice config- 
urations under natural circulation flow with respect to 
leakage reactivity effects in a coolant cell upon sodium 
voiding, mass flow rate and pressure drop in the coolant 
channels, and linear power output. The basis for comparison 
of the three lattices was equal fuel rod size and volumetric 
composition in the cores. The open-hexagonal lattice v;as 
found to be preferable from the standpoint of leakage reac- 
tivity upon sodium voiding because of its larger coolant 
channel nuclear hydraulic radius. For the same cell bulk 
temperature rise the mass flow rates per cell flow area, 
u5 , , /A ,,, were found to be equal for the three lattice 
flow cells. The pressure drop through the lattices were 
approximately equal; the small differences in pressure drop 
were due to the differences in spacer configuration. The 
linear power outputs of the three lattices were found to be 
equal for the three lattices assuming equal core pressure 
drop. Example calculations of the above quantities were 
presented. It was also found that the fluid flov; and heat 
transfer characteristics of the three lattices were equal 
under forced circulation for both laminar and turbulent 
flow. Finally, the applicability of the equivalent diameter 
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of 

conceot for use in calculation heat transfer and fluid flow 
characteristics of the open-hexagonal flow cell was examined. 
It was found that the conventional expression for equivalent 
diameter can be used for the open-hexagonal lattice. 

6.2 Discussion 

The open-hexagonal lattice was found to be superior 
to the triangular and square lattices v/lth respect to leakage 
reactivity upon sodium voiding. If more sophisticated 
calculations confirm that the added negative leakage com- 
ponent of the sodium void reactivity effect in the open- 
hexagonal lattice is indeed several dollars, there would be 
a definite advantage from a safety standpoint in employing 
this lattice in either natural or forced convection LMFBR's. 

However, from the standpoint of hydraulic and thermal 
performance examined in this study, the open-hexagonal lat- 
tice was not found to be superior to the triangular or 
square lattices in any \-iay. The triangular lattice v/ould, 
in fact, be preferable to an open-hexagonal lattice since 
more fuel pins can be packed into a given area of core while 
still maintaining prudent pin to pin spacings. Assuming the 
present FFTF pin spacing of 0.0625 inches for a triangular 

lattice and equal pin diameters of d = 0.300 inches for 

P * 

both lattices, the fuel pins in the open-hexagonal- lattice 
would have to be touching per Eq. (B-1) in order to retain 
an equal fuel volume fraction in the tvjo lattices. If a pin 
to pin spacing for the open-hexagonal lattice were to be 
held at 0.040 inches, the as calculated in Appendix C 
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would have to be held to around Q.306 without considering 

stainless steel swelling and to O.2O8 v;ith consideration 

given to stainless steel swelling. Larger pin sizes for 

the open-hexagonal lattice would allow larger V^, but the 

increase in would not be significant if the pin to pin 

spacing xvas held at a reasonable value 0.020 inches). 

This can be seen from the decreasing slope of the versus 

d curve of Fig. C.l. Current and future design values for 
P 

are found in Table C.l. 

Since the natural convection performance of all three 
lattices is equal in the hydraulic and thermal aspects 
examined, improvements in the linear power output achieved 
by all three lattices can be suggested using Eq. (A- 10 ), 



2 1218 

gpp Z D * n^ ^ At , c„ T ^ 

q- = ( - r)0.555(_)1.555(^)1.0^ 



0.092la.°’®N f, 
e s 



(A-10) 



n. 



In Chapter 4 a value for the ratio of the source/ 

sink elevation difference and the effective primary circuit 
length, of 0.5 was employed. This ratio can be increased. 

Using values from the FFTF design proposals for a 400 Mvjt 

/ m 

LMFBR'' a typical value for is 



N 



e 



AP 

pumpj ^ 

AP 

core 



core 



11.6 ft , 



where 

AP 

pump 



145 psl 
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AP . = 100 psi 

core 



L = 8 ft 

core 



By employing a Z of up to say 2U feet, rhe ratio Z /N 

G G C 

could be made as large as 2.0. For a value of Z /N = 2.0, 

e e 

the linear power output calculated in Chapter 4 would be 
Increased by a factor of (2. 0/0. = 2.l6. In Chapter 
4 the value of the core temperature rise was 400 °P. If 
this value could be raised to 450 °F, the linear power out- 
put would be increased by a factor of (450/400)^*^^^ = 1.2. 

A core At of 400 ^F is about as large a value as is prac- 
tical. A larger value of At may produce a thermal shock 
upon reactor scram v;hich the primary system could not with- 
stand. The linear power output could also be increased by 
spacing the fuel pins farther apart vjhile increasing pin 
diameter. The wider spacing would increase the equivalent 
diameter, D , and coolant flow area, S, thus increasing q' 
since q’ = f ( (D^ * = f(D°*^^^S). The core 
temperature rise. At, would also increase with this change 
in spacing through the relation ’q'riy£, = q - <^c At. Hovjever, 
the spacing contribution probably vjould not increase q' by 
much. In Chapter 4 the core length, 1, v;as four feet. The 
length 1 could possibly be reduced by as much as two feet, 
thus raising q' by the factor (4/2)^**^ = 2. Combining the 
above suggested contributions for increasing q’ gives a 
total factor of about 5.2. With these contributions 'q' 
would equal 4.5^1 kvj/ft instead of 0.8/3 kw/ft as calculated 
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in Chapter 4. However, 4.54 kw/ft is in reality a peak 
value for the natural circulation case. It is less than 
50 ^ of the desired 10 kw/ft average and less than 30/ of 
the 16 kw/ft peak (assuming axial and radial peaking factors 
of about 1.25 and 1.26) needed to match projected 1000 Mwe 
forced convection plant operating conditions ^ Unles.s a 

better combination of the variables discussed above can be 
achieved, full power operation under natural circulation is 
not practical. 

A summary of the parameters assumed and calculated in 
this study is presented in Table VI. 1. 

6.3 Future Work 

Since the leakage reactivity effect in the open-hexagonal 
lattice provides a significant additional negative reactivity 
contribution compared to the triangular and square lattices, 
a m.ore sophisticated calculation by Monte Carlo or other 
m.ethods should be made to determine a precise value for the 
sodium void negative reactivity component in the open- 
hexagonal lattice. The model used in this study yielded 
only approximate answers under the assumption of an infinite 
core. An experim.ent measuring the relative and absolute 
negative reactivity of the lattices should also be consid- 
ered. However, since the reactivity change is on the order 
of one percent of the effective multiplication factor, the 
necessary precision for measuring the reactivity change 
v/ould be difficult to design into the experiment. 
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Models of the three lattices could be constructed to 
test the pressure drop performance of the different config- 
urations. Freon df water could be used as the test coolant. 

6.4 Conclusion 

This study has shown that from hydraulic and thermal 
aspects, the triangular, open-hexagonal and square lattices 
are essentially equivalent under forced and natural circu- 
lation. For natural circulation operation it was shown 
that the linear pov/er output of the three lattices v-jas far 
less than the desired value for full power operation, even 
under an optimum combination of constraints. Therefore, 
nome of the lattices are practical for full power operation 
under natural circulation. Furthermore, unless the open- 
hexagonal lattice is found to contribute significantly to the 
negative component of the sodium void reactivity effect, 
the triangular lattice is preferable for forced convection 
operation. 
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Appendix A 



Linear Povjer Relationship for Natural Circulation 

in a Reactor 



If z is the elevation above a datum surface (see Fig. 
A.l) at which the pressure is constant, the difference in 
pressure between elevations z^ and z^ in a static fluid 
(no friction nor acceleration) is 




(A-l) 



Integrating completely around a closed circuit which has 
the fluid in ti-jo legs at different densities gives 









dz + 



—A. 

p^dz). 
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Introducing the temperature coefficient of volumetric 
expansion, p: 



^h ” ^c ^c^c^^h ” ^c^ * 

gp 3 

APg = -g^ j (th - • <"'-2) 

^ t 

Usj.ng o’ = wc.^(tj^ - tj^) \Mhere q is the heat produced per 
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Fig. A.i, Thermal circulation driving heads. 
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unit cell, and oJ is 'the mass flow rate per cell and 

Zj - = Zg = Zq + (Z^ + Z|^)/2, 



APb = 



gPePeq^e 

ScV 



(A-3) 



From the definition of the Fanning friction factor. 




2H pZJ? 



(A-4) 



an expression for the friction pressure drop, is 

obtained as follows: 



APb 



2 '2 

Sc°eP ® 



(A-5) 



where f , the factor accounting for spacer pressure drop, 
s 

is defined by Eq, (3-10), and u = a>/pS. Equating APg and 
APg and solving for cu gives 



(jj 






(A-6) 



In this derivation is the total effective length of the 
primary coolant system: 



AP 

N = — 



AP 



core 



core 



v;h'ere L is the length of the core. The quantity N is 
core e 

a measure of the equivalent length of the complete primary 

circuit (core, pipes, and heat exchanger) expressed in 

terms of an extended core of height N and channel diameter 

e 

D^, Substituting the relationship fp = 0.046Re”*^*^ 

= 0. 046ijl®* ^/(D^ w/S)^‘ ^ for turbulent flow into Eq. (A-6) 
gives 



UJ = 



crP,n^7 

S q^o.357 

V OP ' 

0*09ai^-Vpfs 



(A-7) 



Using the relationships q = n^q'l, where n^ is the number 

of effective fuel rods per unit cell, and "q = i^c^At, where 

At = t. - t^: 
h c 



q' = 



u»c At 






(A-8) 



Substituting Eq. (A-7) into Eq. (A-8) yields 



_ 0.357 Cp^t 

o f 

^ e p s 



(A-9) 



Solving for q', Eq. (A-9) becomes 



P 1 P 1 8 

gPp Z D S n.^ ^ c-c-i- At , c„ , p, 

i' = ( !1)0-555(_,1.555(,J.)1.0, 



(A-10) 



- 75 - 



Refer to Fig. 3.6, which is reproduced below as Fig. A. 2, 
for the physical illustration of these results. Fig. A. 3 
shows a comparison of the driving forces in the forced and 
natural convection cases. 

For natural convection the following observations are 
made: If it is desired to increase q' , this is done by 

allowing an increase in At. This in effect increases the 
driving head allowing a greater flow rate, cu per Eq. (A-8) 
at which the AP = AP„. The combination of an increase in 

ID K 

both At and u3 increases q according to the relation 
q = and thereby "q’ increases for fixed n^ and 

1, since "q' = q/n^l. 

Note if At is fixed, i. e., a maximum allowable core 
temperature drop is specified, the cell operating condition 
for constant cell parameters n^ and 1 is fixed at the inter- 
section of the three lines At, APj^, and APg. This inter- 
section fixes w and hence "q and “q' through "q = ^fixed^p^^h'^c^ fixed 

and "q' = "q/n 1 = wc At/n 1 from Eq. (A-10). 
r p r 
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a. Intersection points given by Eq, {A-1) , 




b. Eq. (A-IO) yields the ivhlch corresponds to the 
intersection points of the three lines; i. e., for con- 
ditions of AP^ = Ap and At fixed. 

SD K 

Fig. A. 2. Natural circulation operating conditions. 
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a. Relationship between driving pressure and mass 

flow rate in a forced convection system where, for a fixed 

At, = uJc At/n 1, and co depends on the pump pressure, 
p r 




b. Relationships betvjeen driving pressure, mass flow 
rate, and core temperature rise in a natural convection 

sj'stem where, for a fixed At, q' = wc At/n 1 is found by Eq , 

- . 

(A-10), and w is determined as the intersections of lines 
of At, and 



Fig. A. 3. Comparison of driving forces in forced and 
natural- convection. 
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Appendix 3 



Derivation of Hydraulic Diameters as Functions of 
Pit^ and Diameters of Fuel Pins 



Refer to Fig 1.1 for finding fuel volume fractions, V^: 



2ra^A 

W)h = -rr- 



TTd 



6^P^^j3 



3j3p^ 



9 Ph 



(B-1) 



fVM = ? P 

1 p p5 

2P5J3 



Id?/4 TTd 



2 

p _ 



3j3pf 



6 p. 



(B-2) 






TTd^/4 

P 



^ "sq 



(B-3) 



In the derivations of the linear power output capa- 
bilities of the. lattices, the assumptions made were that 
the same fuel pin diameter was to be used in each lattice, 
and the fuel volume fractions v;ere to be equal. Hence, v/e 
obtain the following relations betv/een lattice pitches: 



(nK = (''Pt 



and 



TT 



js.f = 
3J3 Ph 



TT d ^ 

2j3p' 



P?/Ph = 3/2 . 



(B-4) 
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In like manner 



Psq'^Ph = 



(B-5) 



By definition D = 



4 X cross sectional flov; area 
wetted perimeter 



(De)h = 






2TTd 



3j3Ph - % 

7rd 



(B-6) 






'*(|pt(Pt/2)^ - 



2j3p? - TTd^ 



iTTd 

2 p 



ird 



(B-7) 



Using Eq. (B-4), 






2-ra(|pS) - ird^ 



■ 2^h 

ird 



^ = (De>h ' 



(B-8) 



(°e)sq 



4(p^ - TTd^/4) 4p^ - TTd^ 

^ sq p" ’ _ ^sq p 



TTd 



TTd 



(B-9) 



Using Eq. (B-5)/ 






e ' sq 



TTd 



(B-10) 



The above results shoiv that for equal pin diameters and 
fuel volume fractions, the equivalent diameters of the 
three lattices are equal. 



i' 
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Appendix C 
Lattice Parameters 

Table C.l lists some of the lattice parameters v/hlch 
are being used in current large LliFBR's and which are 
proposed for future large LMFBR’s. Table G.2 and Fig. C.l 
show volume fraction of fuel and triangular lattice pin 
pitch and gap spacing based on a gap spacing of 0.040 inches 
for the open-hexagonal lattice (the minimum assumed prac- 
tical) and equal pin diameters for both lattices. 

The equations used in these calculations are Eq. (B-1) 
and Eq. (B-2). The calculation procedure was to select dp 
and, using 6^ = 0.040 Inches, compute p^ and from 

Eq, (B-1). Then taking (Vp^^ = (V^)^, p^ and 5^ are computed 
from Eq. (B-2) for the same assumed pin diameter. Fig. C.l 
shows that as dp decreases, 6^ approaches 6^ with an accom- 
panying drop in V^. For reasonable pin sizes as shown, and 
constraints of constant d and VL between lattices, the 
following conclusions are drawn: 

1. The pin to pin gap for a triangular lattice 

( 0.120 in.) is considerably larger than that for an 
open-hexagonal lattice (0.040 in.). 

2. The fuel volume fraction for the lattices with 
reasonable pin diameters (0.250 to 0.300 Inches) are 
about 0.45 to . 0 . 47 . Significant Increases above these 
values require large (and uneconomical) increases in 
pin diameter. 



Current and Proposed LT-tPBR Parameters 
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B 



Table C.l. continued 




Table C.l. continued 
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Table C.l. continued 
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Table C.2 

Pin to pin spacing and pin pitch for triangular 
and open-hexagonal lattices which have equal 
pin diameters and equal volume fractions of fuel. 



(in.) 






(in-) 


Ph (in.) 


P- (in.) 


0.25 


o.o 4 o 
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0.290 


0.355 


0.30 


0.040 
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1 1 

0.340 
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C 

• 

0 


0.538 


0.45 
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0.490 
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Triangular lattice pin to pin spacing, 6^ (in) Fuel volume fraction 
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Fig. C.l. Assumed and calculated parameters for the open 
hexagonal and triangular lattices. 



Triangular lattice pin pitch, p, (in) 
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3, The fuel volume fraction, V^, used In the above 
calculations was based on a single fuel cell In a fuel 
subassembly and was equal to 1 - sodium volume 

fraction). Taking Into account the volume of the core 
taken up by the clad, the gap between clad and fuel, 
the subassembly housing, the space between subassem- 
blies, and the control assemblies, the actual fuel 
volume fraction, V^, for the whole core Is 0.306, as 
will be shov/n In detail below. 

Refer to Fig. C.2 for calculation of the actual core 

fuel volume fraction, V^. The fuel subassembly contains 

217 fuel pins of 0,300 Inches outside diameter and has a 

of 0.472, Based on a unit length of fuel assembly, the 

Inside volume of the fuel subassem.bly Is 217V„ , . /Vi 

IU0JL pin X 

3 

= 32.2 In. , Dividing the subassembly Into six equilateral 
triangles, the height, h, of a triangle becomes one half 
the Internal distance across the flats. Using a wall thick- 
ness of 0.150 In. and one half the subassembly to subassem- 
bly spacing of 0,025 In,, the subassembly pitch divided 
by tvjo Is h' = h 0,175 in. One equilateral triangle of 
the subassembly has a volume of g-(32.2) = 5.37 In. . The 



volume of one triangle equals ^bh = 2 ^^)^ “ “ 5.37 



In.'^ or h = 3.04 In. Now h' = 3.04 + 0.175 ^ 3.215 In. 

The volume each subassembly occupies, In the core based 

on the distance from subassembly center across the flats 
to the midpoint betv;een subassemblies Is = 6(^b'h') 
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Flg. C.2. Core subassembly configuration 
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35.75 in. . Using a pin diameter 



- 6(~(i^h')h') = 2j3h'^ = 
of 0.300 in., a clad thickness of O.OI 5 in., and a fuel to 
clad radial gap of 0.002 in., the volume of fuel in a sub- 
assembly is = ( 217 ) Q - ! . 9 . 3^1 = 12.06 in.^. 

Assuming there is one control assembly for every ten fuel 
subassemblies, the actual fuel volume fraction in the core. 



12.06 

V. — = ( ) = 0.306 

^SA-^^^SA 35.75 + 3.575 

This number is low compared to the fuel volume frac- 
tions (V^ = 29.2 to 46.7) of the reactors in Table C.l. 

Our value of does not take into consideration the 
irradiation induced swelling of the stainless steel sub- 
assembly walls. Because of this sv/elling the subassembly 
to subassembly spacing must be 0.25 inches instead of 0.05 
inches. Therefore, h’ = 3.04 + 0.275 = 3.315 in., and 

O 

VgA = 38.07 in.-^. Consequently, our value of is lower: 

= 0 . 288 . The fuel volume fractions of the cores in 
Table C.l do not reflect the stainless steel swelling. The 
methods used to calculate the volumetric compositions of 
the cores for the reactors in Table C.l v;ere not explained 
in the references cited. Therefore, it is difficult to make 
a good comparison of the volumetric com.positions of the 
actual cores in Table C.l and our theoretical core. Hov/ever, 
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if the of 0.306 is too low for practical use, an increased 
pin diameter is feasible vjith the use of uranium carbide 
fuel which has a higher melting point than does uranium 
dioxide fuel. Some designs are proposing pin diameters 
of 0.400 inches for future 1000 Mv/e LMFBR's fueled v;ith 
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Appendix D 

Relative Linear Power Output for Various Lattice 
Configurations Under Forced Circulation 

The average linear power output of a lattice flow cell 
is given by 



v;here 



q’ 



uic At 

, 

n„l 



(D-1) 



to = pAV 



Assuming the pumps are moving the coolant at the same total 

rate through the three lattices, and in vievj of the fact 

that the volumetric compositions of the three lattices are 

equal, the coolant velocity will be the same in the three 

lattices. It is shovm in Appendix B that the equivalent 

diameters of the three lattices are equal. The lattices 

have the following numbers of effective rods per flow cell: 

(n ), = 2: (n ). = 1/2: (n ) =1. 

\ r'h •’ ' r't ■’ ' r-'sq 

Since 



D X (v;etted perimeter) 



(D-3) 



'h 



D^(27Td^) Dud^ 
e ^ - G p 
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4 2 

D (ifrd ) D 7rd , 

e^2 p^ _ e P _ 1 a 

— lx “ n 



(d-4) 






(D-5) 



I 
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A 



sq 



D (TTd ) 

4 




(D-6) 



Using the above expressions for flow areas in Eq. (D-2), 
which is substituted into Eq. (D-1), gives the following 
relative values for the linear power output of the three 
lattices : 



^ ^ pA^VCpAt/81 ^ A^/2 ^ ^ 







The above ratios show that the linear power output for the 

three lattices are equal under forced circulation for 

eaual V' and d . 

f p 

The pressure losses through the different lattices are 
the same since 

4lpV^ 

AP = f , (D-7) 

D 2g 
e c 

and all the parameters in Eq. (D-7)^ including fp and are 
equal for the three lattices. The corresponding coolant 
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velocities, V, are shoxvn to be equal in the lattices in 
Tables IV. 1 and IV. 3. The forced convection heat transfer 
coefficients are also the same for the three lattices since 

Nu = hDg/k = + Cg(Pe)'^, 

Nu = + C2(DgVpCpA)'^^ (D-8) 

and all the values in Eq. (D-8) are again equal for the thre 
lattices, including the D^'s as discussed in Appendix B. 

From the above ratios and correlations it is seen that 
the linear poxver output, pressure drop, and heat transfer 
characteristics are the same for the three lattice config- 



urations. 
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Appendix F 
Nomenclature 

(unless defined as they appear) 
Meaning 

cross sectional floxv area 
heat capacity 
diameter 
friction factor 
gravitational constant 
heat transfer coefficient 
thermal conductivity 
multiplication factor 
length 
a quantity 
pipe length 
Nusselt number 

average pov/er output per cell 

average linear povjer output per pin 

pitch 

perimeter 

Peclet number 

pressure drop 

radius 

Reynolds number 

cross sectional flow area 



Units 

2 2 
ft^, in. 

BTU/lb-°P 

ft, in. 

ft/hr^ 

BTU/hr-ft^-^F 

BTU/hr-ft-°P 

ft, in. 

ft 

kw 

kw/ft 

in. 

in. 

Ib/ft^, psi 
in. 

ft^, in. 



2 



i 

i 






Units 



English Meaning 



At 


axial temperature difference per 
(inlet to outlet) 


cell 


V 


average coolant velocity 


ft/hr 




fuel volume fraction = 1 - V., 

Na 

v;here V,. Is the sodium volume 
Na 

fraction 






Actual fuel volume fraction in a 
real core 


— 


Z 


height 


ft 


z 


elevation difference 


ft 


Greek 






P 


fraction of delayed neutrons 




P 


temperature coefficient of 
volumetric expansion 


0 

1 


6 


pin to pin spacing 


in. 


5 


boundary layer thickness 


in. 




dynamic viscosity 


Ib/ft-hr 


V> 


kinematic viscosity 


ft^/hr 


P 


density 


Ib/ft^ 


O' 


ratio of areas, A^/A^ 


— 


CO 


average mass flow rate per flow 


Ib/hr 



cell 



Subscripts 

B due to buoyancy forces 

e equivalent 

em empirical 

entrance, exit 



E 



i 






Subscrl 

g 

h 

hy 

P 

r 

R 

s 

SA 
■ sq 
t 

th 

T 

¥ 

1 



bs Meaning 

grids 

open-hexagonal lattice 

hydraulic 

pin 

rods 

due to friction 
spacers 
subassembly 
square 

triangular lattice 

theoretical 

total 

vjetted 

avjay from gridded region 
in gridded region 
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